Additivity of decoherence measures for multiqubit quantum systems 
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We introduce new measures of decoherence appropriate for evaluation of quantum computing 
designs. Environment-induced deviation of a quantum system's evolution from controlled dynamics 
is quantified by a single numerical measure. This measure is defined as a maximal norm of the 
density matrix deviation. We establish the property of additivity: in the regime of the onset of 
decoherence, the sum of the individual qubit error measures provides an estimate of the error for a 
several-qubit system. This property is illustrated by exact calculations for a spin-boson model. 
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Dynamics of open quantum systems has long been a 
subject of study in diverse fields Recent interest 

in quantum computing has focused attention ^ |^ 
on quantifying environmental effects that cause small 
deviations from the isolated-system quantum dynam- 
ics. During short time intervals of "quantum-gate" func- 
tions, environment-induced relaxation/decoherence ef- 
fects must be kept below a certain threshold in order 
to allow fault-tolerant quantum error correction . The 
reduced density matrix of the quantum system, with the 
environment traced out, is usually evaluated within some 
approximation scheme, e.g., 0. In this work, we intro- 
duce a new, additive measure of the deviation from the 
"ideal" density matrix of an isolated system. For a sin- 
gle two-state system (a qubit) this measure is calculated 
explicitly for the environment modeled as a bath of har- 
monic modes 0, e.g., phonons. Furthermore, we estab- 
lish that for a several-qubit system, the introduced mea- 
sure of decoherence is approximately additive for short 
times. It can be estimated by summing up the deviation 
measures of the constituent qubits, without the need to 
carry out a many-body calculation, similar to the ap- 
proximate additivity expected for relaxation rates of ex- 
ponential approach to equilibrium at large times. 

In most quantum computing proposals, quantum al- 
gorithms are implemented in the following way. Evolu- 
tion of the qubits is governed by a Hamiltonian consist- 
ing of sin gle- qubit operators and of two-qubit interaction 
terms |9|, Il0| . Some parameters of the Hamiltonian can 
be controlled externally to implement the desired algo- 
rithm. During each cycle of the computation the Hamil- 
tonian remains constant. This ideal model does not in- 
clude the influence of the environment on the computa- 
tion, which necessitates quantum error correction. The 
latter involves non-unitary operations |Q] and cannot be 
described as Hamiltonian-governed dynamics. 

The accepted approach to evaluate environmentally in- 
duced decoherence involves a model in which each qubit 
is coupled to a bath of environmental modes 0. The 
reduced density matrix of the system, with the bath 
modes traced out, then describes the time-dependence 



of the system dynamics 0, 0]- Because of the in- 
teraction with environment, after each cycle the state of 
the qubits will be slightly different from the ideal. The 
resulting error accumulates at each cycle, so that large- 
scale quantum computation is not possible without im- 
plementing fault-tolerant error correction schemes [6lll2l|. 
These schemes require the environmentally induced de- 
coherence of the quantum state in one cycle to be below 
some threshold. Its value, defined for uncorrelated single 
qubit error rates, was estimated to be between 10~^ 
and 10-4. 

To evaluate error rate for a given system, one should 
first obtain the evolution of its density operator. In realis- 
tic cases various approximations are utilized A i JjJ ■ 
The most familiar are the Markovian-type approxima- 
tions used to evaluate approach to the thermal state 
at large times. It has been pointed out recently that 
these approximations are not suitable for quantum com- 
puting purposes because they are usually not valid 0, 0| 
at low temperatures and for the short cycle times of quan- 
tum cornputation. Several non-Markovian approaches 
0, 0, 0, [3 have been developed to evaluate the short- 
time dynamics of open quantum systems. 

Additional issues arise when one tries to study deco- 
herence of several-qubit systems. One has to consider the 
degree to which noisy environments of different qubits are 
correlated If all the qubits are effectively immersed 
in the same bath, then there is a way to reduce deco- 
herence for this group of qubits without error correction 
procedures, by encoding the state of one logical qubit in 
a decoherence free subspace of the states of several phys- 
ical qubits [lil Hif . In a large scale quantum computer 
consisting of thousands of qubits, it is more appropriate 
to consider qubits immersed in distinct baths, because 
these errors represent the "worst case scenario" that ne- 
cessitates error-correction. 

Once we have the density matrix p{t) for a single- or 
few-qubit system evaluated in some approximation, we 
have to compare it to the ideal density matrix p^'^ (t) cor- 
responding to quantum algorithm without environmental 
influences. We have to define some measure of decoher- 
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ence, D{t), to compare with the fault-tolerance criteria, 
which are not specific with regards to the error-rate def- 
inition pj. It is convenient to have the measure non- 
negative, and vanishing if and only if the system evolves 
in complete isolation. Since explicit calculations beyond 
one or very few qubits are exceedingly difficult, it would 
be also useful to find a measure which is additive (or at 
least sub- additive), i.e., the measure of decoherence of a 
composite system will be the sum (or not greater than the 
sum) of the measures of decoherence of its subsystems. In 
this work, we identify two measures of decoherence which 
have these desirable properties. We prove the asymptotic 
additivity and derive explicit results for decoherence in a 
short-time approximation. 

Typically, the Hamiltonian of an open quantum sys- 
tem interacting with environment has the form H = 
Hs + Hb + Hi, where Hs is the internal Hamiltonian, 
Hb is the Hamiltonian of the environment (bath), Hj 
is the system-bath interaction term. Over gate-function 
cycles, and between them, the terms in H will be con- 
sidered constant. Therefore, the overall density ma- 
trix R{t) of the system and bath evolves according to 
R{t) = exp {-iHt/h)R{0) e^p (iHt/n). Usually the ini- 
tial density matrix i?(0) = yo(0)®6 is assumed to be a 
direct product of the initial density matrix of the system, 
p(0), and the thermalized density matrix of the bath, O. 
The reduced density matrix of the system is obtained by 
tracing out the bath modes, p{t) = Tib R{t). 

To quantify the effect of the interaction with the 
bath 17], we consider the deviation, cr(t), of the 
reduced density matrix from the ideal p'"^\t) = 
exp {—iHst/h)p{0) exp (iHst/h), 

a(t)=p(t)-p«(t). (1) 

As a numerical measure of the deviation we will use the 
operator norm \\ct\\^. It is defined as the maximal 
absolute value of the eigenvalues of the matrix a. For a 
2x2 deviation matrix, ||cr||;^ = (|crii|^ + |cri2|^)^^^- 

Note that ||cr||;^ is not only a function of time, but it 
also depends on the initial density operator p(0). Due to 
decoherence, it will increase from zero as time goes by. 
However, its time-dependence will also contain oscilla- 
tions at the frequencies of the internal system dynamics, 
as illustrated in Fig. 1. The effect of the bath can be 
better quantified by the maximal deviation norm, 

D{t)= sup {\\a{t,p{0))\\^). (2) 

P(0) 

One can show that < D{t) < 1. This measure of de- 
coherence will typically increase monotonically from zero 
at < = 0, saturating at large times at a value D(oo) < 1. 
The definition of the maximal decoherence measure D{t) 
looks rather complicated for a general multiqubit sys- 
tem. However, in what follows we show that it can be 
evaluated in closed form for short times, appropriate for 
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FIG. 1; Schematic plot of the deviation norm for a two-level 
system interacting with the Ohmic bath P] of bosonic modes 
in the short-time approximation 0. The oscillatory curves 
illustrate the ||o-||_^ norms for three different initial states p(0). 
The monotonic curve shows the D norm, obtained as ||o-||^ 
maximized Q over p(0). Note that in this case D{oo) — 1/2. 

quantum computing, for a single-qubit (two-state) sys- 
tem. We then establish an approximate additivity that 
allows us to estimate D{t) for several-qubit systems as 
well. 

To be specific, let us outline the calculation of this mea- 
sure of decoherence of a two-level system (spin, qubit) 
interacting with a bosonic bath of modes. Thus, we take 

Hs = -ihn/2)a,, HB = Ek^kala^, (3) 
Hi = (TxJ2ki9k'^i+ g*kH)- (4) 

Here ujk are the bath mode frequencies, a^,, aj, are the 
bosonic annihilation and creation operators, gk are the 
coupling constants, ax and ctz are the Pauli matrices, 
and hQ > is the energy gap between the ground (up, 
|1)) and excited (down, |2)) states of the qubit. 

The dynamics of the system can be obtained in closed 
form in the short-time approximation Q . The expression 
for the density operator of a spin-1/2 system is 

P'tnn 

^ gj[(iS<,+-E„--Ep-B„)t/2-B2(t)(,,^-r,„)V4]^ 

Here the Roman-labeled states, with j = m,n,p,q = 
1,2, are the eigenstates oi Hs corresponding to the eigen- 
values Ej/h = {~iyn/2. The Greek-labeled states, |C), 
are the eigenstates of ax, with the eigenvalues 77^ = C, 
where C = = ±1- The spectral function B^{t) is 
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defined according to 

B^t) = 8Efel5fc|'wfe ' 8m\LOkt/2) coth(/3wfc/2). (6) 

Here (3 = h/k^T, Ub is the Boltzmann constant, and T 
is the initial bath temperature. This function has been 
studied in With pi2(0) = |/9i2(0)|e*'^, we get 

IkWII, = i[l-e-^^(*)]{[pn(0)-p22(0)]^ 

+ 4|pi2(O)|^in2[(17/2)t + 0]}'/'. (7) 

In Fig. 1, we show schematically the behavior of ||(T(t)||^ 
for three representative choices of the initial density ma- 
trix /o(0). Generally, ||o'(i)||j;^ increases with time, reflect- 
ing the decoherence of the system. However, oscillations 
at the system's internal frequency are superimposed, 
as seen explicitly in iQ. 

The effect of the bath is better quantified by 
the maximal operator norm D{t). In order to 
maximize over /9(0), we can parameterize p(0) = 
U {P\1){1\ + {1- P)\2){2\)U\ where < P < 1, and 
U is an arbitrary 2x2 unitary matrix, 

/ e*("+T) cos 61 e^'^^'-'iUinO \ . . 

\ -e^(^-")sin0 e-'("+^)cos0 ) ' ^ ' 

By using the rigorous definition 0| of the norm, \\A\\\ — 
sup^_^q((V'|^^^|V')/(V'IV')): one can show that it suffices 
to put P = 1 and search for the maximum over the re- 
maining three real parameters a, 7 and 9. The result, 
sketched in Fig. 1, is 

D(t)^\[l-e-^'(% (9) 

In quantum computing, the error rates can be signif- 
icantly reduced by using several physical qubits to en- 
code each logical qubit [ig. Therefore, even before ac- 
tive quantum error correction is incorporated evalu- 
ation of decoherence of several qubits is an important, 
but formidable task. Consider a system consisting of two 
initially unentangled subsystems 6*1 and S'2, with deco- 
herence norms Ds^ and Ds^, respectively. We denote the 
density matrix of the full system as PS1S2 ^-nd its devia- 
tion as (J 3^821 ^-nd use a similar notation with subscripts 
Si and S2 for the two subsystems. If the evolution of sys- 
tem is governed by the "noninteracting" Hamiltonian of 
the form iJgjSa = + , where the terms Hs^ , Hs2 
act only on variables of the system 6*1, S2, respectively, 
then the overall norm Ds^Sa '^^-'^ be bounded by the sum 
of the norms Ds^ and Ds2 , 

B>S^S2 = SUplla^^sJlA = SUp|lps^5^ - = 
p(0) p(0) 

sup 1 1 ®Ps, - Psl ''Psl I U = sup 1 1 <7sfPs, + Psl ®CTS2 I U 
p(0) p(0) 

< sup\\as^»psj\\ + sup\\p^sl^asj\x 

p(0) p(0) 

< sup ||crsillA+ sup \\crs2h = ^Si+ Ds2- (10) 

PSi(O) PS2(0) 



In general, initially unentangled qubits will remain 
nearly unentangled for short times, because they didn't 
have enough time to interact. Therefore, the inequal- 
ity DsiS2 ~ + Ds2 is expected to provide a good 
approximate estimate for DsiS2i i-^-j the measures for 
individual qubits can be considered approximately addi- 
tive. For large times, the separate measures become of 
order 1, so this bound is not useful. Instead, the rates of 
approach of various quantities to their asymptotic values 
are approximately additive in some cases. In the rest of 
this work, we focus on the short-time regime, of no sig- 
nificant interaction effects, and establish a much stronger 
property: We prove the approximate additivity for ini- 
tially entangled qubits whose dynamics is governed by 
(I3I4II with coupling to independent baths. 

For intermediate calculations, we use the "diamond" 
norm introduced by Kitaev [l^. which we denote by 

K{t) = \\T- r«|U - sup||{[r - r«]«^G}f?l|Tr. (ii) 

e 

Here, the trace norm is given by 11^11^^ — Tr V A. 
The two superoperators are defined as generating the ac- 
tual and "ideal" dynamical evolutions, 

T{t):p{Q)^p{t), T«(t):p(0)K-p«(i), (12) 

Ig is the identity superoperator on a Hilbert space G 
whose dimension is the same as that of the corresponding 
space of the superoperators T and T^^\ and g\s a. density 
operator in this space of twice the number of qubits. 

Consider again the composite system consisting of the 
two subsystems Si, S2, with the "noninteracting" Hamil- 
tonian HsiS2 — +Hs2- The evolution superoperator 
of the system will be Ts^s2 — Tsi^Ts2, and the ideal one 

will be T^'^o = T^'^»T^*^ The diamond measure for the 
0102 01 02 rr^ 

whole system can be expressed [l^ as 

\\iTs.-4>T,^+4'!HTs2-T^;^)\\o 

< II (T,^ - 4>T,J. + \\4>iTs2 - 4'% = 

- T^ifiur^jio + iiT«iuiir^^ - il;'iu = 

\\Ts, - 7^5?IU + \\Ts, ~ T^S2\\o - ^5. + Ks,. (13) 

The step that follows the inequality, in H13|l . proved in 
[1^ as the "stability" property, was established specifi- 
cally with the definition Ijlll) . 

The approximate inequality K 3-1^82 ~ ^Si + ^82 for 
the diamond norm K has thus the same form as for the 
norm D. Let us emphasize that both relations apply 
assuming that for short times the subsystem interactions, 
directly with each other or via their coupling to the bath 
modes, have had no significant effect. However, there 
is an important difference that the relation for D further 
requires the subsystems to be initially unentangled. This 
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restriction does not apply to the relation for K. This 
property is particularly useful for quantum computing, 
which is based on qubit entanglement. However, even in 
the simplest case of the diamond norm of one qubit, the 
calculations are extremely cumbersome. Therefore, the 
measure D is favorable for actual calculations. 

In general, one can prove that 2D{t) < K{t) < 2. For 
a single qubit, calculation within the short-time approx- 
imation actually gives K{t) = 2D{t). Since D is gen- 
erally bounded by K/2, it follows that for the specific 
model considered, with a bosonic heat bath acting inde- 
pendently on each of the qubits, the multiqubit norm D 
is approximately bounded from above by the sum of the 
single-qubit norms even for the initially entangled qubits: 

D<K/2< Y.,K,/2 = E,(l-e-^.^)/2 = , (14) 

where j — 1, . . . ,N labels the qubits. 

For short times, one can also establish a lower bound 
on D{t). Consider a specific initial state with all the 
qubits excited, p(0) = (|2) (2|)i» . . . »(|2) (2|)Ar. Then ac- 
cording to (O, p{t) = /9i(i)» . . . »p7v(i), where pj{t) = 
(l/2){[l-e-^?(*)]|l)(l|+[l+e-^'(*)]|2)(2|}. The right- 
bottom matrix element of the (diagonal) deviation oper- 
ator, cr2«,2"(0 = -1 + 2"^nj[l + e"-^?^*^], can be ex- 
panded as (1/2) Ej + o(Ej ^|(0), because B{t) is 
linear in t for small times. The largest eigenvalue of a{t) 
cannot be smaller than (T2".2"(0- It follows that 

D > hEjB^+oiE.B^) = E,D,+o{J:^D,), (15) 

where we used © for short times. 

By combining the upper and lower bounds, we get the 
final result for short times, 

m^E,D,it)+o{j:^D,it)). (16) 

This completes the proof of additivity of the measure 
D{t) at short times, for spin-qubits interacting with 
bosonic environmental modes. This result should ap- 
ply as a good approximation up to intermediate, inverse- 
system-energy-gap times, e.g., and it is consistent 
with the recent finding '2Q| that at short times decoher- 
ence of a trapped-ion quantum computer scales approxi- 
mately linearly with the number of qubits. 
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